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We investigate light meson bound state contributions to six quark Green's functions and establish 
that the latter has poles corresponding to baryon bound states. We estimate light baryon masses 
including the Roper resonance. Constituent quark model and the Parton model are natural con- 
sequences of this quantum field theoretic picture. We elaborate this analysis to pentaquarks and 
heptaquarks. An almost model independent prediction for the mass of 6 pentaquark is 1.57Geu. 

PACS numbers: 11.55.-m, 11.80.-m, ll.80.Jy, 12.38.-t, 12.39.Mk, 12.40.Yx, 14.20.-c 
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We look for a baryon pole in six quark Green's func- 
tions (qqqqqq) in QCD. Non-perturbatively this has many 
contributions, in particular let us concentrate on contri- 
butions due to quark-quark scattering. In a relativis- 
tic Quantum held theory like QCD this can be inferred 
from quark-antiquark scattering by crossing symmetry. 
We assume the color singlet s channel of qq scattering 
does have color singlet meson poles which are the ob- 
served pion, p,ui,K etc(FigpP). The wave- function <f> n of 
these mesons are realised as eigensolutions of the Bethe- 
Salpeter(BS') equation wherein a kernel is subsumed to 
exist from the underlying theory, QCD. 
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where <f> n obeys the BS equation (FigPJj]). M n is the 
mass of the meson of type n. 



K = S(p) S(p) / Kcj> 



(2) 



S(p) and S(p) are the exact propagators for the quarks. 
The BS wave-function does have a spectral representa- 
tion which follows from causality |l||2||. It is given by 
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d/3 4>n(y,(3) 
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The spin angular momentum, flavour details are sup- 
pressed in the above equation. More precise definitions 
along with spin-orbit coupling details are presented in 
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FIG. 1: (a)The Meson contributions to the quark anti-quark 
scattering kernel, (b) The Bethe-Salpeter equation for the 
Meson 



[2( along with asymptotic behaviour of the BS wave- 
function which are deduced from the structure of the 
BS equation. For completeness we state that there is a 
generic spectral representation for the quark propagator 
satisfying the Schwinger-Dyson(S'D) equation in QCD 
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We have introduced two important quantities|3j fh and 
S n . rh is the threshold mass i.e., only for momenta p 2 > 
m 2 will the quark propagator have an imaginary part. It 
may in general be a pole if the quark is not confined or a 
branch cut if it is confined. Similarly S n is the relativistic 
analogue of the size(in units of mass) of the meson of 
mass M n . In particular the meson wave- function <p n is 
a function of two complex variables, p 2 and p 2 . Only 
for momenta p 2 > S 2 or p 2 > 5„ the wave-function <p n 
has singularities. The threshold mass m is in principle 
determined form the SD equation while the meson size 
S n and meson mass M n are self-consistently determined 
from the BS equations. From eq.JSJ it is easy to see that 



< S n < min(rh, rh) 



(5) 



The lower inequality follows from the the fact that all 
amplitudes have no singularities in the space-like do- 
main(here it amounts to saying the meson has a finite 
size) while the upper bound is self-evident from the BS 
equation, eq.© (fh and fh being the threshold masses of 
the quark and the anti-quark). A more precise estimate 
of S n demands the knowledge of the exact structure of 
the kernel in the BS equation. From the definition of 
<5„(eq.(|3J)), we note that a bound quark in a meson has a 
threshold singularity at 5 n which is less than that of the 
unbound quark(eq.j5jl). 

Consider the amplitude of three quarks going to three 
quarks. We notice that there are contributions due to 
three different quark-quark scattering. One such contri- 
bution of interest is the one shown in Fig[2] ■ This involves 
one internal loop momentum k as explicitly shown in the 
Figj2j. By choice we take qi and rj(with r\ — 0) as inde- 
pendent momenta while the external momenta p L and pi 
follow from momentum conservation at the vertices. To 
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FIG. 2: Baryon pole(constituent quark picture). 



avoid the ambiguity of double counting the quark prop- 
agator, we introduce the kernel used in the BS equa- 
tion. Using spectral representation for the BS wavefunc- 
tion(eq.©) this diagram is completely well defined and 
the loop momentum, k, integration can be done. We are 
interested in the singularity structure of this graph. We 
apply the Hadamard/Landau conditions Q on all vari- 
ables being integrated over, which includes yi and ft 
coming from each of the spectral representations of the 
meson vertices. To make the computation transparent 
first notice that we have to consider terms like 

A. = Vi(k + nf - (1 - yi )pi - w (l - Vi)M* - Sf - ft (6) 

This is equivalent to an effective mass term for the quark 
propagator which has momentum (fc + r^). Landau equa- 
tions demand that ft = and ^ = 0, 1 or ^ = 0. 
Furthermore looking for the lowest effective mass with re- 
spect to p 2 i.e., Tjp- = we require yi — 1. Consequently, 

we find A » (k + r^) 2 — S 2 . To evaluate the nature of the 
singularity it is convenient to employ Cutkosky cuts[5j as 
shown by the dashed line in Fig|2]. With suitable condi- 
tions on cj)(y, (3) it yields six Dirac-delta functions which 
impose the following constraints on the momenta 

= ^Mf + ^-kf (7) 
(r° + k ) = ] /5? + (r l + k)2 (8) 

Four of the Dirac-delta functions will eliminate the loop 
momentum k integration completely One will be the 
overall pole condition on the total energy(eq.^H) and 
the other Dirac-delta function restricts the external mo- 
menta. The loop momentum k integration also yields a 
self-consistency condition(Landau equation) namely the 
vanishing of the following four momenta 

£ (-aKft- AO + &<(»■< + *;)) =0 (9) 

i 

for some a,,6j > 0. The physical interpretation of this 
equation is that if the masses of all lines in the graph of 
Fig|2] are taken to be real, then the process depicted in 
this graph can occur in real space-time classically pro- 



vided they obey this constraint jg. The total energy is 



e = 53a? + r?) =53 v^ 2 + (Qi -k) 2 + \Af + w + fc) 2 

»=1 i=l 

(10) 

Minimising E, all spatial momenta vanish and we get the 
lowest baryon mass Mg, 

3 

Mfl = 53(M, + *i) (11) 

i=l 

For this case, the conditon in eq.|0 can be satisfied triv- 
ially. The mass of the lightest baryon, the proton, is in 
terms of the mass of the pion taken to be M n = 139Mev 
and the size of the pion, S w , inferred to be YJAMev{\fm). 

The picture that we have deduced from this analysis 
is precisely that of the constituent quark model. The 
interpretation of the graph(Fig[2|) is that there is a pole 
corresponding to a baryon state \B) and corresponds to 
the following matrix element, 

(qqq\B)——± T —(B\qqq) (12) 

S IVl — (— Id 

where s — E 2 in the rest frame of the Baryon. Conse- 
quently we infer the overlap function (B\qqq) in terms of 
the explicit meson wave-funcitons. Furthermore for the 
lightest baryons the three quark momenta pi are such 
that they are constrained to be pf = (313Mei>) 2 . In the 
rest frame of the baryon these energies of the quarks de- 
compose into masses of the mesons and 5 W (bound 
quark mass in the meson) as shown by the cut in Fig[2]. 

In the theory Fig[2] is not the only pole contribution 
to the six quark Green's function. Indeed there are many 
others. Firstly this graph can repeat several times by it- 
self or repeat with other quark propagations in between. 
All these also have singular Dirac-delta function contri- 
butions corresponding to the baryon pole. We envisage 
a systematic expansion wherein these contributions re- 
sult in a finite renormalisation of the baryon mass Mb- 
Therefore our estimate for S n — HAMev subsumes these 
contributions. In 't Hooft's -h expansion[7j where N is 
the number of colors, we know that every additional me- 
son function implies a factor of -4= . Some of these finite 

renormalisations are controlled in the j? expansion. In- 
deed these contributions makes the evaluation of masses 
of excited baryons tedious and also detailed knowledge 
of meson wave-functions becomes necessary. 

In addition there are other corrections which allows 
us to infer Feynman's parton model. For example, us- 
ing the same technique we can look at Green's functions 
where the three quark state goes to three quarks and sev- 
eral other gluons and quark anti-quark pairs as shown in 
FigEl- All these graphs also have the same pole contribu- 
tion as shown by the cut in Fig[3] and indeed this graph 
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FIG. 3: Baryon pole(Parton picture), 
corresponds to 

(qqq\B) L. (B\qqq, gg(etc)) (13) 

tS IVl Q ~\~ 1£ 

where g stands for gluons. From this we infer 
(B\qqq, gg(etc)) matrix elements explicitly. That is the 
baryon state has a finite probability amplitude to pro- 
duce various other degrees of freedom in additon to the 
valence quarks. 

(B\ = (B\qqq)(qqq\ + (B\qqqgg) {qqqgg\ + ■■■ (14) 

From our analysis we see that each term on the the r.h.s 
of the above equation has the same mass eigenvalue Mp 
as dictated form Fig[3]. It should be noticed that the 
constitutent quark picture of Fig[2] is disturbed because 
these constituents can also exchange energy as shown in 
the Fig[2| leading to parton distribution functions. We 
also note in the conventional parton picture, quarks with 
momenta Pi(Fig[2]) emit gluons but we find in general 
there are quarks with lower energy, of the order of Si 
which also emit gluons as shown in Fig Although 
formally all these are precisely given in terms of meson 
wave- functions and quark and gluon propagators, a sys- 
tematic convergent expansion is an open problem. 

Stability of the baryon state would require 
that Mb < 3m. We have no a priori estimate 
of m in QCD but a self-consistent picture to be pre- 
sented later shows that m < 591Mev. In a model theory 
such as crQCD 8] it has been estimated^ for u,d, s 
quarks to be about 550Mev. 

Excited baryons can come about because of excited 
mesons in the intermediate states of Fig [2]. Note that 
excited mesons such as p, lu etc are necessarily unsta- 
ble becuase they do decay into 2tt or 3tt states. Con- 
sequently the meson pole is in the complex energy do- 
main. Our analysis of the Green's functions can be car- 
ried out even in the complex plane consistent with causal- 
ity (i.e., with Feynman it prescription) namely posi- 
tive energy pole has negative imaginary part and vice 
versa. The external momenta pi and pi then have to 
be complex. As a consequence the width of the baryon 
is necessarily greater than the widths of these interme- 
diary mesons. For example if we consider intermediate 
mesons to be p and n with masses M p and 2M n , we 



get M B = 2(M 7r + 5 n ) + [M p + 5 p ) » 1474Mew for 
d p « ^ (since the size of the p meson has to be larger 
than that of the the pion) . This matches with the Roper 
resonance and we would also find that the widths obey 
r Roper > ~ 150Mei>, which is consistent with phe- 
nomenology. As we see in this example predicting higher 
states demands knowing meson properties, M n and 5 n . 
The latter being unavailable we have to develop theoret- 
ical models. As another example, consider the A particle 
which is made up of u, d, s quarks. The intermediate 
mesons in Fig[2] have to be two it and K°. From the 
masses of these mesons we get the mass of the A baryon, 
M A = 1.2Gev for S K w 4*. 

Our method of estimating the spectrum of baryons al- 
though unconventional is consistent with quantum field 
theory. It is also consistent with the Hamiltonian eigen- 
value problem. The complication is that the eigenvalue 
problem is stated in the basis of infinite Fock states. It 
is interesting to note that the same eigenvalue problem 
fea. (|14f> ) can be solved by using Green's functions and 
their singularities. We can explicitly get the mass eigen- 
values and baryon states in terms of meson states. Now, 
it has to noted that if we take a system which is essen- 
tially non-relativistic, i.e., the threshold masses are large 
compared to the binding energy, bi , then the mass of the 
2-body meson bound state is given by, Mi ks 2m — bi 
and Si ~ am where a is some interaction coupling con- 
stant. Consequently we will estimate the baryon mass 
from Fig|2] to be 3(2m + am — bi) which is much greater 
than 3m (mass of three free quark states) implying that 
the system is unstable. This shows our procedure is in- 
applicable for all heavy quark systems. Similar consid- 
erations on charm and u, d quark systems also results in 
instability for the baryon bound state. 

Exact chiral symmetry is another limit where we can 
draw conclusions. This as discussed in Q is a singular 
limit. Due to Goldstone's theorem we have an exact re- 
lationship between the pion and the quark propagator 
namely 4>tt{p) oc ^TrS(p) showing that d v = m and 
also Mn — 0. Consequently the mass of the baryon is 
Mb = 3m implying that the baryon is almost unsta- 
ble. Hence in the constituent quark model estimate, the 
baryon cannot be bound. However there are finite renor- 
malisation effects alluded to earlier. These effects can 
give rise to a weakly bound baryon state. 

Now consider the recently discovered pentaquark 
candidate which has been discussed in the literature [llT|. 
This is a ududs color singlet system with mass 1540Mev. 
In our picture the graph and cut shown in Fig^] gives the 
lowest pentaquark possible. The mass of the pentaquark 
9 is therefore 

Me =3(M n +S„)+M K +M 7T = Mp+Mk+M* = IhllMev 

(15) 

where Mp is the mass of the proton or neutron and Mk 
is the mass of the K meson. Interestingly this prediction 
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FIG. 4: Pentaquark pole. 



does not explicitly depend on 6^ if the mass of the proton 
is given. Figpl] shows that the spin of this state is \ as in 
the case of the proton. As for as the width is concerned, 
since it is made up of only stable mesons, it is expected 
to be narrow. 

We may look at other pentaquarks such as ududu, our 
prediction for its mass will be Mp + 2M„ = VlVTMev 
and spin i, which is close to the mass of the A baryon. 
This leads us to suspect that in the standard A reso- 
nance which has a width of about YlQMev there exists 
a pentaquark resonance as well. The standard A reso- 
nance is the uuu state which in our picture comes about 
from symmetry considerations in the same way as in the 
constituent quark model. In terms of graphs, there are 
actually two graphs of the Fig[2] type and the cancel- 
lations between them yields a higher mass to the uuu 
system than an uud system. Similarly the masses of hep- 
taquark states such as udusdsu can also be predicted to 
be M P + 2M K + 2AU = 2313Meu and this resonance will 
be a little broader since it has many decay channels. 

In our considerations here, the important infrared 
properties of the quarks and mesons that are relevant to 
the description of baryons are the threshold masses, m, 
meson masses, M n and sizes, S n . We have experimental 
knowledge of M n and to a lesser extent of 8 n inferred from 
form factors. In principle fh and S n can be inferred from 
Lattice simulation of QCD. Another attempt is to model 
QCD qq scattering kernel in the BS equation(FigP3])in 
a reliable manner. crQCD 8] is one such approach where 
the kernel is taken to be purely a one gluon exchange 
and a singular propagator. Here q is the exchanged 
momentum between the quark and the anti-quark. This 
was analysed in the large N expansion with g 2 N(g is 
the gauge coupling constant) also small. In this theory 
which is asymptotically free, we find spontaneous sym- 
metry breakdown of chiral symmetry and PCAC0, El- 
The threshold mass fh of u, d, s quarks is numerically es- 
timated to be 550Mev and 5 n can be estimated. 

In the above discussion even if the quark is not confined 



i.e., the quark propagator has a pole at fh, as long as the 
inequalities alluded to in the paper are satisfied the re- 
sults remain true. Color non-singlets such as diquark sys- 
tems can also be investigated in our picture. The graph 
of this system is that of Fig[2] but with with the middle 
quark line removed. In 4 dimensions the loop momentum 
is uniquely fixed and there is no residual Dirac-delta func- 
tion to restrict the total energy of the Green's function. 
This implies the diquark system has no pole but only a 
branch cut singularity. This may have important con- 
sequences in scattering processes. Similarly if we take a 
four quark system, they can be bound like the three quark 
system with a mass ^(m v + S n ) even though they can 
never be a color singlet. Looking at non-singlets one can 
consider a qqq system which has the same color quantum 
numbers as a fundametal quark. This yields the same 
graph as in Fig0] but without two quarks. The estimate 
of the mass of this state being Zm^ + S w = 591Mev and 
is necessarily a branch cut. Consequently consistency of 
our picture demands that the threshold mass be, 

m < 591Mev (16) 

Color singlet criterion has played very little role in our 
analysis, if at all these considerations have to come from 
other theoretical aspects of QCD. It is assumed that 
mesons are color singlets. It should be borne in mind that 
production of these multi-quark states becomes harder 
because these have to obey the criteria of eq.© in the 
real phase space. 
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